Structure of multiphoton quantum optics. I. Canonical 
formalism and homodyne squeezed states 
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We introduce a formalism of nonlinear canonical transformations for general systems of multi- 
photon quantum optics. For single-mode systems the transformations depend on a tunable free 
parameter, the homodyne local oscillator angle; for n-mode systems they depend on n heterodyne 
mixing angles. The canonical formalism realizes nontrivial mixings of pairs of conjugate quadra- 
tures of the electromagnetic field in terms of homodyne variables for single-mode systems; and in 
terms of heterodyne variables for multimode systems. In the first instance the transformations yield 
nonquadratic model Hamiltonians of degenerate multiphoton processes and define a class of non 
Gaussian, nonclassical multiphoton states that exhibit properties of coherence and squeezing. We 
show that such homodyne multiphoton squeezed states are generated by unitary operators with 
a nonlinear time-evolution that realizes the homodyne mixing of a pair of conjugate quadratures. 
Tuning of the local-oscillator angle allows to vary at will the statistical properties of such states. 
We discuss the relevance of the formalism for the study of degenerate (up-)down-conversion pro- 
cesses. In a companion paper, "Structure of multiphoton quantum optics. II. Bipartite 
systems, physical processes, and heterodyne squeezed states", we provide the extension of 
the nonlinear canonical formalism to multimode systems, we introduce the associated heterodyne 
multiphoton squeezed states, and we discuss their possible experimental realization. 

PACS numbers: 03.65.-w, 42.50.-p, 42.50.Dv, 42.50.Ar 



I. INTRODUCTION 

The study of the nonclassical states of light has re- 
cently attracted renewed attention because of the key role 
they may play, beyond the traditional realm of quantum 
optics, in research fields of great current interest, such 
as laser pulsed atoms and molecules pj; Bose-Einstein 
condensation and atom lasers |2( ; and quantum informa- 
tion theory Fundamental physical properties for an 
efficient functioning of the quantum world such as inter- 
ferometric visibility, robustness of superpositions against 
environmental perturbations, and degree of entanglement 
are typically enhanced by exploiting states which exhibit 
strong nonclassical features. 

The simplest archetypal examples of nonclassical states 
are of course number states, whose experimental realiza- 
tion is however difficult to achieve. Moreover they share 
very few of the coherence properties that would be desir- 
able both in practical implementations and in fundamen- 
tal experiments. The most important experimentally ac- 
cessible nonclassical states are the two-photon squeezed 
states HHH]. The y 

are Gaussian, exhibit several im- 
portant coherence properties, and can be obtained by 
suitably generalizing the notion of coherent states y|. 
Squeezed states can be easily introduced by linear Bo- 
goliubov canonical transformations and the associated 
eigenvalue equations for the transformed field operators 



b(a,rf): 6|* >p= /?|* > . They can be produced in 
the laboratory through the dynamical evolution of para- 
metric amplifiers y|, and provide a useful tool in vari- 
ous areas of research. For instance, they have been pro- 
posed to improve optical communications @ and to mea- 
sure and detect weak forces and signals like gravitational 
waves |lC|. Moreover, twin beams of bipartite systems, 
i.e. two-mode squeezed states, are maximally entangled 
states, a property of key importance in quantum compu- 
tation and in quantum information processing. 

Realistic and scalable schemes of quantum devices and 
operations with continuous variables might however re- 
quire the realization of multiphoton processes. In this 
respect, it is crucial to investigate the existence and 
structure of multiphoton nonclassical states of light. A 
challenging goal is to define suitable multiphoton gen- 
eralizations of the effective Hamiltonian description of 
two-photon down-conversion processes and two-photon 
squeezed states. Natural candidates should be nonclas- 
sical states obtained by nonlinear unitary evolutions as- 
sociated to anharmonic Hamiltonians and multiphoton 
down-conversion processes. In turn, nonlinear unitary 
evolutions might be of great importance in the imple- 
mentation of universal quantum computation with con- 
tinuous variables (CV) systems 0,0]- Experimental 
realizations could be obtained by considering the dynam- 
ics of the polarizability in nonlinear optical media |l3[ : 
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where P is the polarization vector, x the n-th order 
susceptibility tensor, and E the electric field. Implemen- 
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tation of higher order multiphoton parametric processes 
involves several terms in the expansion Eq. For in- 
stance, fc-photon parametric down conversions involve 
all contributions at least up to the term with coupling 
^( fc ) j whose strength in nonlinear crystals is in general 
extremely weak for k > 2. It must however be remarked 
that coherent atomic effects, such as electromagnetically 
induced transparency and coherent population trapping 
manipulations of photons in cavities provide new promis- 
ing techni ques to generate large and lossless optical non- 
linearities 

Phcnomcnological theories of multiphoton paramet- 
ric amplification, based on expansion Eq. Q, were in- 
troduced by Braunstein, Caves and McLachlan [l5j . 
who considered nonlinear interaction terms of the form 
Zkd' — zIcl producing fc-photon correlations. The prob- 
lem was numerically addressed by the authors, who 
showed that these interactions generate squeezing and 
display remarkable phase-space properties. Another in- 
teresting approach, but rather abstract since it involves 
infinite powers of the canonical creation and annihilation 
operators, was put forward in ref. (16j where generalized 
multiboson, non Gaussian squeezed states were intro- 
duced. A crucial question left unanswered by the above- 
mentioned attempts is that, although the first nonlinear 
order in expansion Eq. (JIJ can be associated, through the 
linear Bogoliubov transformation, to an exactly diagonal- 
izable two-photon Hamiltonian and to exact two-photon 
squeezed states 0, 0, higher order nonlinearities have 
not been associated so far to exact Hamiltonian mod- 
els of multiphoton effective interactions in a simple and 
physically transparent way. 

In a series of recent papers ^3 113 1 a first step was re- 
alized in this direction by defining canonical transforma- 
tions that allow the exact diagonalization of a restricted 
class of multiphoton Hamiltonians. In this formalism one 
adds to the linear Bogoliubov transformation a nonlinear 
function of a generic field quadrature. The canonical con- 
ditions impose very stringent relations on the parameters 
of the transformations, and the resulting Hamiltonians 
describe a very peculiar and restricted class of nonlinear 
interactions, not easily amenable to realistic experimen- 
tal realizations. 

Stud ying the simplest case of a quadratic nonlincar- 
ity [TtI Il8|. one determines multiphoton squeezed states 
both in the case of nonlinear functions of the first quadra- 
ture X\ = (a + a^/Vz and in the case of nonlinear func- 
tions of the second quadrature Xi — —i(a — a))/y2. 
These states exhibit interesting nonclassical statistics 
and squeezing in the quadrature associated to the non- 
linearity; they may be denoted as single- quadrature mul- 
tiphoton squeezed states (SQMPSS). The unitary opera- 
tors associated to the two transformations are a compo- 
sition of squeez ing , displacement, and a nonlinear phase 
transformation 01 (see also The scheme, although 

limited to one-mode systems, still provides some insights 
for experiments in quantum information exploiting mul- 
tiphoton processes |2(j. In fact, the single-quadrature 



multiphoton squeezed states include as a particular case 
the generalized "cubic phase" states (originally intro- 
duced in the framework of quantum computation [TT| ^1 
proposed by Bartlett and Sanders by adding displace- 
ment and squeezing to the pure cubic phase transfor- 
mation |12| . The single-quadrature canonical formalism 
(and the associated single-mode, single-quadrature mul- 
tiphoton squeezed states) is thus very limited because it 
amounts only to a pure (nonlinear) phase transforma- 
tions on a single quadrature. Moreover, it does not allow 
nontrivial extensions to multimode systems and nondc- 
generate processes. 

In the present and in a companion paper, that we shall 
denote as Part I and Part II, we show that these difficul- 
ties can be overcome and that it is indeed possible to 
introduce a general canonical formalism of multiphoton 
quantum optics. In the present paper (Part I) we de- 
termine the most general nonlinear canonical structure 
for single-mode systems by introducing canonical trans- 
formations that depend on generic nonlinear functions of 
homodyne combinations of pairs of canonically conjugate 
quadratures. The homodyne canonical formalism defines 
a class of single-mode, homodyne multiphoton squeezed 
states; it includes the single-mode, single-quadrature 
multiphoton squeezed states as a particular case; and 
introduces a tunable free parameter, a local-oscillator 
mixing angle, which allows to interpolate between differ- 
ent multiphoton model Hamiltonians and to arbitrarily 
vary the field statistics of the states. In the companion 
paper "Structure of multiphoton quantum optics. 
II. Bipartite systems, physical processes, and het- 
erodyne squeezed states" (Part II) |2lJ we extend the 
multiphoton canonical formalism to multimode systems. 
We show that such extension is realized by nonlinear 
canonical transformations of heterodyne combinations of 
field quadratures. The scheme defines a structure of mul- 
timode, heterodyne multiphoton squeezed states that re- 
duce to the homodyne states for single-mode systems. 
In Part II we also show that the heterodyne squeezed 
states and the associated effective multiphoton Hamil- 
tonians can be realized by relatively simple schemes of 
multiphoton conversion processes |2l|. In this way we 
introduce a complete hierarchy of canonical multipho- 
ton squeezed states: a) multimode heterodyne squeezed 
states; b) single-mode homodyne squeezed states; and c) 
single-mode, single-quadrature squeezed states. In the 
limit of vanishing nonlinearity the multiphoton squeezed 
states reduce to the standard (single-mode or multimode) 
two-photon squeezed states. 

To construct a general canonical formalism of multi- 
photon quantum optics one must first circumvent the re- 
strictions following from the canonical conditions. In par- 
ticular, the prescription that a general, nonlinear mode 
transformation be canonical prevents the possibility of 
introducing arbitrary nonlinear functions depending si- 
multaneously on two conjugate quadratures Xi and A 2 . 
In fact, if the form of the nonlinear function is not con- 
strained at all, the canonical conditions force the nonlin- 
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ear coupling to be trivially zero. 

It is however possible to define a general canonical 
scheme by introducing a simultaneous nonlinearity in 
two conjugate quadratures if the nonlinear part of the 
transformation is an arbitrary function of the homo- 
dyne combination v2|?7| (cos#Xi + sin 6X2 ) of the two 
quadratures (77 = |?7|expi0 arbitrary complex number). 
Such a canonical structure allows naturally for a local 
oscillator angle 9 mixing the quadratures. The mixing 
is a physical process that can be easily realized, e.g. by 
a beam splitter positioned in front of a nonlinear crys- 
tal. Tuning the continuous parameter 9 allows then to 
vary the physical properties, and in particular the statis- 
tical properties, of the associated homodyne multiphoton 
squeezed states. 

The plan of the paper is as follows. In Section [n] 
we develop the general formalism of nonlinear canoni- 
cal transformations for homodyne variables. In Section 
IIIII we study the multiphoton Hamiltonians associated to 
the nonlinear transformations, specializing to the case of 
quadratic nonlinearity. We compute the wave functions 
of coherent states associated to the transformations, as 
functions of the mixing angle 9. In Section llVI we deter- 
mine the explicit form of the unitary operators associated 
to the canonical transformations. They are composed by 
the product of a squeezing, a displacement, and a mix- 
ing operator with nonquadratic exponent which combines 
conjugate quadratures. In Section^] we study the statis- 
tical properties of the homodyne multiphoton squeezed 
states. We compute the uncertainty products, and deter- 
mine the condition for " quasi- minimum" uncertainty. We 
then study the quasi-probability distributions, and the 
photon statistics. We show that these properties depend 
strongly on the local oscillator angle 9. In Section lVll we 
summarize our results and discuss the outlook and ex- 
tensions that are developed in the companion paper Part 
II. 



II. NONLINEAR CANONICAL 
TRANSFORMATIONS FOR HOMODYNE 
VARIABLES 

We could naively imagine that the problem of intro- 
ducing nonlinear canonical transformations for a single 
bosonic mode a should be solved by adding to the stan- 
dard Bogoliubov linear transformation an arbitrary (suf- 
ficiently regular) Hermitian nonlinear function F(a,a') 
of the fundamental mode variables a, a^: 



b = /ia + va) + 7-F(a, a)) . 



(2) 



Requiring that the transformed mode b be bosonic, i.e. 
that [b, b^} — 1, and exploiting the well known formulae 

[a,G(a,a^)} = dG(x, y)/dy\ x = a , y=a i , 



[a t ! G(a,a t )] = -8G(x, y)/dx\ x=a ^ =a i , 



then the condition for the transformation Eq. © to be 
canonical reads 



ImI 2 - M 2 + M 2 ^ 1 ] + 
.apt 



(3) 
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da 



, * dF 
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dF 
daJ 



= 1 . 



It is very difficult to determine the most general expres- 
sion of the nonlinear function F, which allows to satisfy 
the condition Eq. J3J). Nevertheless, if we assume that 
the nonlinear function be hermitian, then canonical gen- 
eralizations of the Bogoliubov transformation do exist, 
and the most general expression is in terms of arbitrary 
Hermitian, nonlinear, analytic functions F of homodyne 
linear combinations of the fundamental mode variables. 
It reads: 



b = fia + va* + jF(ifa + rja*) 



(4) 



with rj = \rj\e %e a complex number. Exploiting the func- 
tional dependence of F on the homodyne combination of 
the modes a and a\ one finds that the general relation 
Eq. reduces to the following algebraic constraints on 
the complex coefficients of the transformation: 

ImI 2 -M 2 = i, 

i?e[e ie (/i 7 * -jAy)] =0. (5) 

With the parametrization 

/i = coshr , v = sinhre 1 ^ , 7 = |7|e i5 , (6) 

we can express the canonical conditions Eqs. @ in the 
form of the transcendental equation 



coshrcos(0 — S) — sinhrcos(5 + 9 — <fi) — 



(7) 



Eq. J7J can be solved numerically For instance, given 
some fixed r, 9 and <f> we can find numerical solutions 
for the phase variable 6, which can be used as an ad- 
justable parameter to assure the canonical structure of 
the transformation. Alternatively, we can look for par- 
ticular analytical solutions of Eq. (0): letting (j> = 0, we 
obtain the simplified expression 



tan 9 tan 5 



-2r 



(8) 



For given values of the local oscillator angle 9 this is 
a relation between the phase S of the nonlinearity and 
the strength r of the squeezing; e.g., fixing 9 = ±7r/4, 
we get tan (5 = =Fe _2r . Setting 9 = —S implies instead 
tan S = er r . Of course, Eq. J7J admits infinite solutions 
which correspond to a great variety of nonlinear canon- 
ical operators. We can however select more stringent 
conditions, imposing 



7T 

±-±fc7T 



7T 

±-±/wr. 
2 



(9) 
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with k, h arbitrary integers. This choice allows to satisfy 
Eq. at the price of eliminating one degree of freedom 
from the problem. In conditions Eqs. © it is obviously 
sufficient to consider k = h = 0. From Eqs. JSJ and (J7J 
it is evident that the modulus of rj is irrelevant in the de- 
termination of the canonical constraints of the transfor- 
mations. Therefore, from now on we set \t)\ = In 
this way, the homodyne character of the transformation 
scheme becomes fully evident. We can in fact express the 
transformation Eq. J3J in terms of the rotated homodyne 
quadratures Xg,Pg defined as 



X„ = 



-iS 



/V2 = Xi cos 8 + X 2 sin 8 



Pg = X e+7T / 2 = -Xi sin 8 + X 2 cos ( 



(10) 



with [Xe,Pe] = i. The transformed mode b can then be 
expressed in terms of the rotated mode ag = ae~ lS , or of 
the rotated quadrature Xg, as 



(11) 



b = fiag + va\ + jF(Xg) , 

with the rotated parameters 

pi = ne w ; v = ve~ l9 . (12) 

From Eqs. (JTDJ , (|11|) . (|12[) the canonical conditions 
Eqs. JSJ) can be expressed in terms of the rotated pa- 
rameters as 



IA| 2 -|-| 2 



1 



Re[(prf* - u*j)} = 0. 



(13) 



The form Eq. Ijlll) of the canonical transformation Eq. Q 
allows the straightforward determination of the coher- 
ent states of the transformed modes, and of the unitary 
operators associated to the transformations, as we will 
show in the following Sections. Obviously, when 7 = 
one recovers the standard linear Bogoliubov transforma- 
tions and the structure of standard two-photon squeezed 
states. 



III. MULTIPHOTON HAMILTONIANS AND 
MULTIPHOTON SQUEEZED STATES 

We now consider the Hamiltonians that can be asso- 
ciated to the canonical transformations Eq. Q. We will 
restrict ourselves to the case of quadratic Hamiltonians 
diagonal in the transformed modes: 

H = tfb, (14) 

whose general expression in terms of a, er reads 

H = \v\ 2 + (\^\ 2 + \v\ 2 )a)a+{^u*a 2 + h.c) 



In terms of the rotated quadratures Xg, Pg, and exploit- 
ing the canonical constraints Eq. Q, the Hamiltonian 
Eq. iJTSJ reads 

H = -l + l^+^Xl + ^-^Pg 2 



+ \!\ 2 F 2 {Xg) + -^=7*(/i - v){Pg, F(Xg)} , (16) 

where {•,•} denotes the anticommutator. The Hamilto- 
nian Eq. I |lti|l is further simplified by exploiting the con- 
ditions Eqs. ©: 

H = ^X 2 -l+ e -^-{Pg± V2\ 7 \e ±r F(X e )} 2 . (17) 

Eq. I|17|l is of the same form obtained in Refs. |l7|- 
jl8j| . but with the fundamental difference that the non- 
linearity is now placed on the homodyne-mixed, rotated 
quadrature Xg rather than on a single one of the orig- 
inal X-s. Moreover the mixing depends on a tunable 
free parameter, the local oscillator angle 8. Eq. 117|) 
shows that the variable Xg is squeezed and that its con- 
jugate variable, in the sense of being antisqueezed of a 
corresponding amount, is the "generalized momentum" 
Pg ± V2\j\e ±r F(Xg). The associated quasi-probability 
distributions are then squeezed along a rotated axis, as 
will be shown in the following Sections. 



A. The case of quadratic nonlinearity 

In studying the statistical properties of the coherent 
states associated to the Hamiltonians Eqs. H15f) - l|16|l we 
will specialize to the case of the lowest possible nonlin- 
earity in powers of the homodyne rotated quadratures, 
i.e. we will consider the quadratic form 



F(Xe) = X 2 g 



(18) 



Inserting Eq. I|18fl in Eqs. (|15fl - lltty . the corresponding 
four-photon Hamiltonian reads 

H Ap = A + (A x a) + A 2 a) 2 + A 3 a i3 + A 4 a t4 + h.c.) 



B a^a + Sia t2 a 2 



Ca t2 a + .Da t3 a + h.c. 



where the coefficients Aj, Bi, C, D are 



(19) 



+ (fj,*^F + v*jaF+ l^F^F + h.c. 



(15) 



A = \v\ 2 + || 7 | 2 , A x = i/x* 7 + + e 2 *V 7 , 
M = ^| 7 | 2 + H*v , A, = ^Vl + ^ 2 *V , 
A 4 = ie 4 ^|7| 2 , B = H 2 + IH 2 + 3|7| 2 , 
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Bi = ^| 7 | 2 , C = \e 2i \m* + + M*7 + ^7* , 



D = e 2i6 \y 2 



(20) 



The Hamiltonian Eq. i|19|) describes one-, two-, three- 
and four-photon processes with effective linear and non- 
linear photon-photon interactions associated to degen- 
erate parametric down conversion processes in nonlinear 
media. We see that the Hamiltonian coefficients Eq. (|20() 
crucially depend on the homodyne angle 9, allowing for a 
great freedom in searching for physical implementations 
of multiphoton states by processes associated to higher 
order nonlinear susceptibilities. 



B. Homodyne multiphoton squeezed states 

We define the coherent states associated to the 

Hamiltonian Eq. (|15|l as the eigenstates (with complex 
eigenvalue (3 = |/3|e^) of the transformed annihilation 
operator b: 



b\*)f, = 



(21) 



Choosing the representation ^ p{xg) = (xg\^f)/3 in which 
the homodyne rotated quadrature Xg is diagonal, the 
eigenvalue equation Eq. 12 1H reads 

d xe ^p{xg) = - — , [(/2 + u) xg 

+ V2jF(xg) - V20]V (x e ) , (22) 

where we have used Pg = —id Xg . The general solution of 
Eq. 1221) is 



Wp(x B ) = /Vexp 
where 

N = 



CXg — b 



dyF(y) 



(23) 



Re[a] 



-1/4 



exp 



{Re[c]f 



2Re\a] 



is the normalization, and the coefficients read 
jjL + v v/2~7 y/2(3 



b = 



jl — v 



(1 — V 



fl — V 



(24) 



It can be easily verified that it is always Re[a) > and 
Re[b] — 0. The wave function can then be expressed in 
the form 



~Re[a}~ 


1/4 




exp 


7T 





Rela] 



xg 



exp 



% Im[b] / dyF(y) + 



Im[c 



Re[c] 
Re[a] 



Im[c]xg 



(25) 



By recalling the definition of the parameters Eqs. © and 
the canonical conditions Eqs. ©, Eq. (|25|l reduces to 



^@{xg) = A exp 



= ±2r 



Xg - 



exp 



^^/2e ±r 



| sin(£ - 6)x e - M / dyF{y) 



( 26 ) 

where A = ir 1 / 4 e ±r / 2 . The wave function Eq. I|2t)|) has 
a Gaussian density in the variable xg, with a squeezed 
variance, while both the phase and the density depend 
crucially on the homodyning, i.e. on the local oscillator 
angle 9. We name therefore the states Eqs. <|21[1 - H26[) 
Homodyne Multiphoton Squeezed States (HOMPSS). 

The dependence on the homodyning 9 has important 
physical implications, especially with regard to the sta- 
tistical properties of the HOMPSS. In general the states 
Eqs. i|25|) - (|26[) in the Xg-diagonal representation dis- 
play non Gaussian terms in the phase, such as a cubic 
phase term in Xg in the case of a quadratic nonlinear- 
ity. The nontrivial structure of the HOMPSS emerges 
clearly when writing in terms of the original quadra- 
tures, for instance in the Xi-diagonal representation 
^ p(x) = (x\ty)f3. Specializing to the quadratic nonlin- 
earity F(Xg) = Xg and for a generic angle 9 ^ 0, ir, one 
has 



^p(x) = A/" exp {iax 2 + bx} Ai 



cx 



c 2/3 



(27) 



where Ai[-] denotes the Airy function, and the coefficients 
a,b,c,d are given by 



1 



■ cot ( 



2 

\72~7 sin 3 9 



b = 



jJL — V 



d = 



2V27sin 2 ' 
(A* - P 



87 2 sin 9 7 sin 



The general non Gaussian character of the HOMPSS is 
thus apparent when writing them in the original field- 
quadrature representation. 



C. Reduction to single— quadrature multiphoton 
squeezed states 

When considering the special cases 9 = and 9 = it/2 
the HOMPSS reduce to the Single-Quadrature Mul- 
tiphoton Squeezed States (SQMPSS) previously intro- 
duced in Refs. [l7l ITsj . In these two special cases the 
canonical transformations Eq. reduce to 



b = i±a + va) + jF(Xi) , i = 1, 2 



(28) 



where i = 1 for 9 = and i = 2 for 9 = ir/2. 

The associated coherent states are defined as the eigen- 
states of b with eigenvalue j3\ in the case of zero phase dif- 
ference between [i and i>, and parameterizing (3 in terms 
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of the coherent amplitude a, j3 = fia+va* (a = ai+ia 2 ), 
they read [l7llis|: 



V7 



-2/', 



where 



exp{i[ciXi + e r, 7 J G(a; 4 )]} , i = 1, 2 , 



n = r , r 2 = -r , 

7i = Jm(7) ,72 = ~Re(j) , 



(29) 



G(z) 



ii 



4 0) = V2a, , = -V2a 2 , 
ci = %/2a 2 , C2 = — v2o;i . 

The expression Eq. for the SQMPSS shows that the 
quadrature associated to the nonlinearity is squeezed, 
while, in the chosen representation, the nonlinear func- 
tion F enters in the phase of the wave packet. Explicit 
non Gaussian densities are again realized if we adopt 
the "coordinate" representation (in which Xi is diago- 
nal) when the nonlinearity is placed on X 2 , or viceversa. 



IV. UNITARY OPERATORS 

The expression Eq. H26fl allows to identify, in terms 
of the homodyne rotated quadratures, the unitary op- 
erator Uhom associated to the canonical transformation 
Eq. (nj, such that the HOMPSS are obtained by apply- 
ing Uhom to the vacuum: 



|*>/3 = U hom \0) 



(30) 



where |0) is the vacuum state of the original field mode 
a: a\0) = 0. The unitary operator Uhom then reads 



Uhom = Ug(Xg)Dg(ag)Sg(C,g) , 



(31) 



where 



De{ae) = exp (aga\ - a, 



gag 



is the standard Glauber displacement operator with ag — 

p,*ft - i>(3*, and 



2 9 



is the standard squeezing operator with = re 1 ^ 2e \ 
Finally, the "mixing" operator Ue reads 



U e (X e ) = exp 



ilm[b] / dYF(Y) 



(32) 



where, exploiting the canonical conditions Eqs. J^J, 
Im[b] = ±\/2|7|e ±r , and the abstract operatorial integra- 
tion acquires a precise operational meaning in any cho- 
sen representation. If we specify to the case of quadratic 
nonlinearity F(Xg) — X|, we have: 



U e {X e ) 



exp 



-ilm[b]^- x 2- 3 ' 2 {ae 



= exp 



+ SaWe-* + 3o ta ae M + 3ae~ ,;e + Sa+e 49 )] . 

(33) 

The structure of Ug is elucidated by expressing it in terms 
of the field quadratures: 

Ug(Xg) = exp [-ilmlb}^ 1 (Xf cos 3 8 + Xl sin 3 6 
+ iXfX 2 cos 2 6 sin 9 + 3XiX 2 sin 2 6 cos 9 
+ 3iXiCos 2 6»sin6» + 3iX 2 sin 2 9 cos 9))] . 

(34) 

The unitary operator Ug Eq. I|34|l depends on all powers 
of the conjugate field quadratures up to n + 1 if the non- 
linearity F is a power of order n of the homodyne quadra- 
ture Xg. Ug is a mixing operator: it mixes nontrivially 
the original quadratures, depending on the values of the 
local oscillator angle, except for the special cases 9 = 
and 9 = n/2 treated in Refs. [l7l Hq . where the mix- 
ing disappears and the nonlinearity is a simple power of 
a single field quadrature (either X\ or X 2 ). Obviously, 
the quadratic nonlinearity is only one of a large class of 
possible choices allowed for F in Eq. I|32|l . so that many 
complex nonlinear unitary homodyne mixing of the orig- 
inal field quadratures can be realized. 

In the special cases 6 = and 9 = tt/2 the HOMPSS 
reduce to the SQMPSS. Then, from the choice of rep- 
resentation of Eq. I|29l) one determines the form of the 
unitary operators that produce these particular multi- 
photon squeezed states: 



Ui = ea;p[ie n 7 i G(X i )]-D(a)5(r), 



(35) 



where S(r) is the one-mode squeezing operator for 
( j and D(a) is the Glauber displacement operator. Wc 
see that in these particular cases the nonlinear part of 
the transformation adds to squeezing and displacement 
a pure nonlinear phase term in one of the quadrature. 
We also notice that with the choice F = X\ we recover 
the cubic phase states proposed in Ref. ^3 as a possible 
tool for the realization of quantum logical gates. 



V. STATISTICAL PROPERTIES AND 
HOMODYNE ANGLE TUNING 

A. Uncertainty products 

When considering the statistical properties of the 
HOMPSS we first study the behavior of the uncertainties 
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in the home-dyne quadratures Xg, Pg. Let us express the 
generalized variables in terms of the transformed mode 
operators b and b^: 



1 



X e = -j=[U?-v*)b+(ji-v)tf] , 



(36) 



P e = -=[{fi+i>)tf ~(fl* + t>*)b-2ilm{jl 1 * -v*j)F(Xg)} . 
v2 

(37) 

The above lead to the following expressions for the un- 
certainties: 



A 2 Xg 

A 2 Pg 



v\ 2 



-\Ji + i>Y + 2Im 2 [jl* 1 - £7*] x 

((F 2 )f3 - (F)}) - 2/m[A*7 - *7*] x 
Im[(ji + v)([F,b%] , 



(38) 



where (-)p denotes the expectation value in the HOMPSS 
l^)^, and [•, •] denotes the commutator. It is evident 
that the nonlinearity affects only APg, as the second of 
Eqs. (|38|) explicitly depends from the form of the function 
F. Considering the quadratic form for F and assuming 
the canonical conditions Eqs. ©, Eqs. I|38|) become : 



A 2 X e 
A 2 Pg 



~ 2 



= 2 C 



±2r 



+ e T2r |7| 2 {l + 4|/3| 2 + 4i?e[e- 2ie /3 2 ]} . (39) 



If we now consider the uncertainty product, and if we 
define (3 = |^|e l «, Eqs. £2) yield 

A 2 X e A 2 P g = I + I|7| 2 e^ { i +4 |^| 2 + 4|/3| 2 cos2(^0)} . 

(40) 

It is to be remarked that this last relation attains its 
minimum for £ — 8 = ±-|: 



A 2 XgA 2 Pg 



(41) 



Eq. I|41|) can be seen as a "quasi-minimum" uncertainty 
relation; in fact, although the second term is not exactly 
zero, for small nonlinearities it will be surely very small 
with respect to the first term (i.e. the Heisenberg mini- 
mum), due both to I7I < 1 and to the decreasing contri- 
bution of the exponential for a suitable choice of the sign 
of r. 



B. Average photon number 

We now turn to the calculation of the average number 
of photons (n) = (0*0) in a HOMPSS |*)^. We special- 
ize to the case of a quadratic nonlinearity. We will show 
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FIG. 1: Mean photon number (n) as function of I7I, for a 
HOMPSS with magnitude of squeezing r — 0.8, coherent 
squeezed amplitude j3 — 3, and different mixing angles: 6 — 
(full line); = ? (dashed line); and 6 — ? (dotted line). 



that (n) is strongly affected both by the strength I7I of 
the nonlinearity and by the mixing angle 8. In Fig. ^ we 
study the behavior of (n) as a function of I7I for fixed 
values of the squeezed coherent amplitude [3 and of the 
magnitude r of the squeezing, and for three different val- 
ues of the mixing angle 8. Due to the canonical condition 
<p = 28, there cannot be pure squeezing for 8 7^ even in 
absence of the nonlinearity, and thus at 7 = we have 
different initial average numbers of photons depending 
on the value of 8. The analytic expression for (n) for a 
homodyne transformation with I7I = and (f) = 28 reads: 



| 2 cosh2r - Re[/3* 2 e 2t9 } sinh2r + sinh 2 r . (42) 



We see from Fig. ^ that (n) needs not show a mono- 
tonic behavior as a function of the nonlinearity. It is in 
fact very sensitive to the mixing angle 8, and although 
it eventually always grows for sufficiently large values of 
I7I, it can however initially decrease, depending on the 
mixing angle 8, and then increasing again very slowly 
at larger values of |7|. This behavior suggests that (n) 
will show even more remarkable properties when varying 
8 for different, fixed values of |7|. In Fig. [3 we ana- 
lyze the behavior of (n) as function of 8, at fixed |7|. 
For I7I = we recover the oscillatory behavior of (n) as 
a function of the squeezing phase cj> = 26 of the stan- 
dard two-photon squeezed state. When letting 8 vary at 
fixed finite values of |7|, the oscillations of (n) become 
faster, and both peaks and bottoms quickly rise to very 
large values. Therefore the nonlinearity plays the role 
of a "quantum coherent pump" allowing for very large 
average numbers of photons even in the case of lowest, 
quadratic nonlinearity. 
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C. Quasi— probability distributions and 
phase— space analysis 

We now turn to the study of the statistics of direct, 
heterodyne, and homodyne detection for the homodync 
multiphoton squeezed states \^}p, showing in particu- 
lar how the statistics is significantly modified by the 
tuning of the mixing angle 9. We will consider the 
HOMPSS Eq. (|26|l . corresponding to the canonical con- 
ditions Eqs. 10, and we will specialize to the lowest non- 
linear function F(xe) = Xg. 

We first consider the quasi-probability distributions 
associated to HOMPSS for some values of 9 and com- 
pare them to the corresponding distributions associated 
to the standard two-photon squeezed states. This will al- 
low a better understanding of the behavior of the photon 
number distributions and of the normalized correlation 
functions that will be computed later. The Q-function 

Q(a) = -|<a|*>/3| 2 . (43) 

gives the statistics of heterodyne detection, and corre- 
sponds to a measure of two orthogonal quadrature com- 
ponents (| a) being the coherent state associated to the 
coherent amplitude a = a\ + ia-i). Figs. E] and El 
show three-dimensional plots of the Q— function of the 
HOMPSS, for fixed values of r, (3 and for two dif- 
ferent values of 9. For 9 = x/2, which corresponds to 
the case F = X$, i.e. no mixing, the plot resembles the 
Q— function for a squeezed state, but we can observe a de- 
formation of the basis, curved along the Re[a] = u\ axis. 
For 9 = 7r/3, a case of true mixing of the field quadra- 
tures, the deformation becomes much more evident and 
the function is strongly rotated and elongated with re- 
spect to the Q-function of the two-photon squeezed state. 
Homodyne detection measures a quadrature component 
X\ = ^m(a-\ + a !)> wnere A is a phase determined by the 
phase of the local oscillator. Homodyne statistics cor- 



FIG. 3: Plot of the Q-function, with r = 0.8, hi = 0.4, 
(3 = 3, 9 = for the canonical conditions S — 9 = —J, 

s + e-4> = f. 




FIG. 4: Q-function, with r = 0.8, \j\ = 0.4, (3 = 3, 9 = f , 
for the canonical conditions S — 9 = —^,8 + 9 — 4>= ^. 

respond to project the Wigner quasiprobability distribu- 
tion onto a x\ axis. With the identification A = 9 we plot 
the Wigner quasi-probability distribution for orthogonal 
quadrature components xe and pe: 

W(x e ,pg) = ^ J dye- 2ip °y**p(xe-y)*0(xe+y) . (44) 

In Figs. El and El we show respectively a global projec- 
tion and an orthogonal section of the Wigner function 
for 9 = vr/2, i.e. no mixing, fixed canonical constraint, 
and intermediate value I7I = 0.4 of the nonlinearity. We 
see from Figs. El and El that the Wigner function dis- 
plays interference fringes and negative values, exhibiting 



FIG. 5: W(xQ,p e ), with r = 0.8, | 7 | = 0.4, (3 = 3, 9 = § , for 
the canonical constraint 5 — 6 = —%,5 + 6 — cj>=%. 




FIG. 6: W(xg,p e ), with r = 0.8, \j\ = 0.4, /3 = 3, 9 = § , for 
the canonical conditions 5 — 9 = —^,8 + 9 — cj>=^. 



FIG. 7: Plot of the Wigner function W(x g ,pg), with r = 0.8, 
I7I = 0.4, (3 = 3, 9 = ~, for the canonical conditions 5 — 8 = 
— f, <5 + # — = 



rise, for suitable values of the parameters, to deformed 
intersection areas with the circular crowns associated in 
phase space to number states. In turn, this will lead to 
modified behaviors of the photon number distribution, 
including possible enhanced or subdued oscillations 
as the mixing angle 9 is varied. Moreover, we expect 
that also the second- and fourth- order normalized cor- 
relation functions will strongly depend on 9, showing a 
deeper nonclassical behavior, for instance antibunching, 
in correspondence of values of 9 associated to stronger 
nonclassical features of the Wigner function, like nega- 
tive values and interference fringes. 



D. Photon statistics 



a strong nonclassical behavior. In Fig. we show the 
Wigner function, for the same values of r, |— y | , /3, but 
for 9 = vr/3, a value that realizes a true mixing of the 
field quadratures. We see that the distribution becomes 
strongly rotated and elongated, with a pattern of inter- 
ference fringes and negative values, providing evidence of 
the complex statistical structure of the HOMPSS when 
a true homodyne mixing is realized. The behaviors of 
the quasi-probability distributions suggest the following 
considerations. We recall that the HOMPSS, although 
being states of minimum uncertainty in the transformed 
("dressed") modes (&, 6^), are not minimum uncertainty 
states for the original quadratures Xi and X%; in fact, 
Eqs. I|39[) - (|41[) show that a further term due to the un- 
avoidable statistical correlations adds to the pure vacuum 
fluctuations. This fact is reflected in the rotation and in 
the deformation of the distributions. But we expect that 
this features will affect also the behavior of other sta- 
tistical properties, like the photon number distribution. 
In particular, shape distortions of the quasi-probability 
distributions strongly modify the original ellipse associ- 
ated to the standard two-photon squeezed states, giving 



We begin by analyzing the probability for counting 
n photons, the so-called photon number distribution 
(PND) P(n) = \(n\^)i3\ 2 , in direct detection, and ne- 
glecting detection losses: 



P(n) 



dx 9 (n\xe}(xe\^)f3 
1 



2™n! 7 r 1 / 2 



dx e e 2 H n (xg)^> p{xe) 



.(45) 



Due to the nonlinear nature of the function F, it is in 
general impossible to write a closed analytic expression 
for P(n), which can however be easily determined nu- 
merically. In Fig. El we plot the PND of the HOMPSS 
for various intermediate values of the local oscillator 
9, together with the PND of the standard two-photon 
squeezed states. As foreseen from the previous phase- 
space analysis, the behavior of the PND strongly de- 
pends on the value of the mixing angle 9, and for suitable 
choices of 9, r and I7I, the PND shows larger oscillations 
with respect to the standard two-photon squeezed states 
and the HOMPSS with 9 = 0. Moreover, the oscillation 
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FIG. 8: Pin) for the HOMPSS, corresponding to the canoni- 
cal constraints 8 — 9 — — £ 8 + 9 — <f> — — ~, for different values 
of the parameters: /3 = 3, squeezing magnitude r = 0.8, and 
strength of the nonlinearity | — y ] = (solid line); j3 = 3, r = 0.8, 
I7] = 0.4, = (dotted line); /3 = 3, r = 0.8, | 7 | = 0.5, 9 = f 
(dashed line); = 3^ = 0.5, |-y| = 0.5, 6» = f (dot-dashed 
line) . 



peaks persist for growing n and are shifted for differ- 
ent values of I7I; this behavior is due to the different 
terms entering in the unitary operator Ue(Xg) Eq. (|34"jl . 
which mix in a peculiar way the quadrature operators for 
6>^0,7r/2. 

Concerning the correlation functions, it is interesting 
to study the behavior of the normalized second order 
correlation function 



(«^) 2 ( a t afl )2 



(46) 



and of the normalized fourth order correlation function 
9 (4) (0) = i4^|> (47) 



for the HOMPSS, and determine the different physical 
regimes. In Figs . I§1 and 1 1 01 we compare g^ 2 \0) as a func- 
tion of the squeezing parameter r for the two-photon 
squeezed state and for the HOMPSS at different values 
of 9. We see that also the correlation function g^(0) 
shows the strong nonclassical features of the HOMPSS; 
both the nonlinearity strength |"y| and angle 9 strongly 
influence the behavior of the curves. In fact, the curves 
deviate from the standard form and saturate at lower 
values with respect to TCS. The most significant fea- 
ture is however obtained plotting g( 2 )(0) as a function 
of 9, Fig. ^] Fig. the plots clearly demonstrate 
that it is possible to pass from a subpoissonian to super- 
poissonian statistics. So, tuning 9, we can have photon 
bunching or antibunching and we can select the preferred 
statistics. Similar considerations can be made about the 
fourth order correlation function <? 4 (0), as it can be seen 
by Fig. 1131 where it is shown the dependence from the 
squeezing parameter r, and by Fig. 1141 where instead 



FIG. 9: The correlation function as function of r for the 
HOMPSS, corresponding to the canonical constraints 8 — 9 — 
— ?, 8 + 9 — <f> — — ^, for several choices of the parameters: 
(3 = 3, 7 = (solid line); (3 = 3, |-y I = 0.4, 9 = (dashed line); 



3, hy| = 0.05, 9 
f (dotted line). 



f (dot-dashed line); ,9 = 3, \^\ = 0.5, 




FIG. 10: The correlation function <j' 2 ' as function of r for the 
HOMPSS, corresponding to the canonical conditions 8 — 9 = 
— ?, 8 + 9 — 4> = ■§ , versus the two-photon squeezed states 
(full line), for different choices of the parameters: 1) f3 = 3, 
[7I = 0.05, = ¥ (dashed line); 2) = 3, hi = 0.2, 9 = f 



(dot-dashed line); 3) f3 



n 



0.5, 9 = f (dotted line). 



g 4 (0) is plotted as a function of the angle 9. Also the nor- 
malized fourth-order correlation function shows a partic- 
ular shape for intermediate value of 9 and varying satura- 
tion levels, depending on the parameters of the nonlinear 
term. Moreover, we have four photon bunching or anti- 
bunching in correspondence of different values of 9. 



VI. CONCLUSIONS AND OUTLOOK 

We have introduced single-mode nonlinear canonical 
transformations, which represent a general and simple ex- 
tension of the linear Bogoliubov transformations. They 
are realized by adding a largely arbitrary nonlinear func- 
tion of the homodyne quadratures Xg,Pg. We have in- 
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FIG. 11: The correlation function g^ as function of for the 
HOMPSS, corresponding to the canonical conditions S — 9 = 
§ , 5 + - 4> = f , for: 1) /? = 3, r = 0.5, | -y I = 0.4 (full line); 
2) (5 = 3, r = 0.4, | 7 | = 0.1 (dashed line); 3) /? = 3, r = 0.1, 
|7| = 0.1 (dotted line). 

sp»»3 



2.B 




FIG. 12: The correlation function g^ 2 ' as function of for the 
HOMPSS, corresponding to the canonical conditions 8 — 9 = 
-§, 5 + 61-0 = §, for: 1) /3 = 3, r = 0.8, |i| = 0.1 (full line); 
2) /? = 3, r = 0.5, |7| = 0.4 (dashed line). 



troduced the Homodyne Multiphoton Squeezed States 
(HOMPSS) defined as the eigenstates of the transformed 
annihilation operator. The HOMPSS are in general non 
Gaussian, higly nonclassical states which retain many 
properties of the standard coherent and squeezed states; 
in particular, they constitute an overcomplete basis in 
Hilbert space. On the other hand, many of their statis- 
tical properties can differ crucially from the ones of the 
Gaussian states. In particular, we have shown the strong 
dependence of the photon statistics on the local oscillator 
angle 6. Among other remarkable features, there are the 
possibilities of exploiting the homodyne angle as a tuner 
to select sub-poissonian or super-poissonian statistics, 
and as catalyzer enhancing the average photon number 
in a state. 

The single-mode multiphoton canonical formalism se- 
lects a large number of non Gaussian, nonclassical states, 
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FIG. 13: The correlation function g w as function of r for the 
HOMPSS, corresponding to the canonical conditions S — 9 = 
— f, 5 + 6 — <fi = — f, versus the two-photon squeezed states 
(full line), for several choices of the parameters: 1) (3 = 3, 
|7| = 0.4, = (dashed line); 2) /3 = 3, \y\ = 0.4, = £ 
(dot-dashed line); 3) /3 = 3, | 1 = 0.4, = f (dotted line). 
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FIG. 14: The correlation function g ™ as function of for the 
HOMPSS, corresponding to the canonical conditions S — 9 — 
-f , 8 + 9 - 4> = -f , for: 1) f3 = 3, r = 0.5, | 7 | = 0.1 (full 
line); 2) /? = 3, r = 0.1, \y\ = 0.2 (dashed line). 



including the single-mode cubic phase state, which gen- 
eralize the degenerate, Gaussian squeezed states. On the 
other hand, both from the point of view of modern ap- 
plications, as e.g. quantum computation, and for ex- 
perimental implementations, generalizations to two and 
many modes are of great interest. 

In the following companion paper, "Structure of mul- 
tiphoton quantum optics. II. Bipartite systems, physi- 
cal processes, and heterodyne squeezed states" (Part II), 
we will extend the canonical scheme developed in the 
present paper (Part I) to study multiphoton processes 
and multiphoton squeezed states for systems of two cor- 
related modes of the electromagnetic field. This exten- 
sion is important and desirable in view of the modern 
developments in the theory of quantum entanglement 
and quantum information. In particular, we will show 
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how to define two-mode nonlinear canonical transforma- 
tions and we will determine the associated "heterodyne 
multiphoton squeezed states" (HEMPSS). In the context 
of macroscopic (quantum) electrodynamics in nonlinear 



media, we will moreover discuss the kinds of multiphoton 
processes that can allow the experimental realizability of 
the HEMPSS and of the effective interactions associated 
to the two-mode nonlinear canonical formalism. 
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